arXiv:1506.04841v2 [nucl-th] 8 Sep 2015 


Isoscalar and neutron modes in the El spectra of Ni isotopes and the relevance of 

shell effects and the continuum 

P. Papakonstantinou* 

Rare Isotope Science Project, Institute for Basic Science, Daejeon 305-811, Korea 

H. Hergert 

National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing , MI j882f-1321, USA 

R. Roth 

Institut fur Kemphysik, Technische Universitdt Darmstadt, 6j283 Darmstadt, Germany 

(Dated: September 9, 2015) 

We study theoretically the electric dipole transitions of even Ni isotopes at low energies, using 
the self-consistent quasi-particle random-phase approximation with the D1S Gogny interaction and 
a continuum-RPA model with the SLy4 Skyrme force. We analyze isoscalar states, isovector states, 
and the dipole polarizability. We define a reference value for the polarizability, to remove a trivial 
dependence on the mass number. We compare our results with data and other calculations, with a 
focus on collective states, shell effects, and threshold transitions. 

Our results support the presence of a strong isoscalar transition, with little or moderate El 
strength, as a universal feature of ordinary nuclei. In moderately neutron-rich Ni isotopes, namely 
68 Ni and neighboring isotopes, this transition is found bimodal due to couplings with surface neu¬ 
trons. An adequate treatment of the continuum states appears essential for describing suprathresh- 
old El strength, especially beyond 68 Ni. Very exotic isotopes (N > 50) are found highly polar¬ 
izable, with practically all their El strength in the continuum. The dipole polarizability and the 
neutron-skin thickness are influenced by shell structure in different ways, so that they can appear 
anticorrelated. A comparison with existing results for lighter (Ca) and heavier (Sn) nuclei suggests 
that the so-called pygmy dipole strength is influenced strongly by shell effects and that, partly for 
that reason, its isospin structure depends on the mass region. 

PACS numbers: 24.30.Gd; 21.60.Jz; 21.30.Fe; 21.65.Cd; 


I. INTRODUCTION 

The electric dipole transitions of nuclei at low ener¬ 
gies, below the giant dipole resonance, have been under 
intensifying investigation for more than two decades [1 
3]. Knowing or predicting correctly the concentration 
of dipole strength, which cannot be attributed to the 
low-energy tail of the giant dipole resonance, is espe¬ 
cially consequential in describing nucleosynthesis pro¬ 
cesses [4]. Through its correlations with other observ¬ 
ables, the amount of strength may be used to con¬ 
strain the isovector properties of the nuclear equation 
of state [5, 6]. 

In recent years, the use of hadronic probes and the 
investigation of unstable nuclei have offered new insights 
into the nature and isospin structure of low-energy dipole 
transitions [1]. The extraction of complete spectra on 
both sides of the particle emission threshold also consti¬ 
tutes important progress [7, 8]. The effort is bound to 
continue, both theoretically and experimentally. New re¬ 
sults on exotic nuclei are under analysis and facilities with 
improved capabilities are planned. Moreover, the already 
available dipole spectra are so complex, as functions of 
the excitation energy, mass, and isospin asymmetry, that 
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a unified theoretical interpretation of them has been elu¬ 
sive. New data will be needed to resolve the situation, as 
well as the use of diverse nuclear models. 

Various factors and mechanisms have been proposed 
or are recognised as relevant for generating low-energy 
dipole transitions in medium-mass and heavy nuclei. A 
variety of microscopic models have been employed in re¬ 
lated studies, mainly the quasi-particle or plain random- 
phase approximation (QRPA or RPA) and models con¬ 
sidering phonon coupling. The presence of strength has 
been attributed to the concentration of single-particle 
transitions in the lhu energy region [3, 9-12], an alpha- 
cluster mode [13, 14], a surface-dipole vibration, akin to 
the other known quadrupole and octupole vibrations [15- 
19] and of toroidal structure [20], or, in the case of 
neutron (proton)-rich nuclei, an exotic neutron(proton)- 
skin oscillation [2]. The fragmentation patterns and 
precise strength can be explained by the coupling of 
the dipole transitions to collective phonons [21-24] or 
collisional damping [25]. Is has been demonstrated, 
in studies along isotopic chains, that also shell struc¬ 
ture and deformation influence the low-energy dipole 
strength [17, 19, 26, 27]. We may add that threshold 
particle resonances are a general feature of potential bar¬ 
rier scattering and may also contribute to the photoneu¬ 
tron cross section just above the emission threshold. The 
question now is, whether certain mechanisms dominate 
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in the different mass and asymmetry domains of the nu¬ 
clear chart, and at what energy in relation to the emission 
threshold. 

The answer to the first part of the question may well 
be different in the case of the isoscalar (IS) or isovector 
(IV) channel. That is the conclusion we have reached 
in Refs. [15-19] for the dipole strength below particle- 
emission threshold and in moderately neutron-rich nu¬ 
clei (below critical shell closures). The IS strength is 
attributed to a collective mode of vibration, which is ex¬ 
cited in N = Z nuclei too; the IV strength, largely to 
single-particle transitions. Our conclusions were based 
on QRPA calculations with the Gogny D1S interaction 
and comparisons with experimental data as well as with 
other calculations, along the Ca and Sn isotopic chains 
and in N = Z nuclei. A genuine neutron-skin oscillation 
may be excited either above threshold or in extremely 
neutron-rich nuclei. We note that relativistic models fa¬ 
vor a neutron-skin oscillation in moderately neutron-rich 
nuclei [28], while IS coherence is generally the outcome 
of non-relativistic studies [10]. Such qualitative disagree¬ 
ments can be attributed to the softness of the symmetry 
energy, which is typically predicted to be different by the 
two types of models. 

Pinpointing the mechanisms which control the dipole 
strength and establishing correlations with other observ¬ 
ables often requires extreme care in interpreting theo¬ 
retical results, and of course data. A case in point 
is the correlation between the dipole polarizability and 
the neutron-skin thickness, which was found model- 
dependent [29]. One question which we will address in 
this work is whether the shell structure affects the po¬ 
larizability and the skin thickness in different ways. As 
regards data, it will be very interesting to confirm exper¬ 
imentally the dramatic role of shell closure revealed in 
several theoretical studies. 

The second part of the question we have posed above is 
related to scattering in the particle continuum. In most 
theoretical studies so far the particle continuum is at 
best discretized (for recent exceptions see, e.g., [30, 31]). 
Broad resonances in the vicinity of the particle threshold 
may not be adequately described. This is another issue 
we will address here. 

In the present work we focus on the Ni isotopic chain, 
which represents an intermediate mass region, between 
the lighter Ca and heavier Sn nuclei. Known isotopes 
span the whole pfg $/2 neutron shell. They include the 
long-lived, symmetric 56 Ni isotope and several stable 
ones. Photoresponse data exist for the latter, below par¬ 
ticle threshold, and for the neutron-rich 68 Ni, above par¬ 
ticle threshold, while more are expected to be analysed 
soon. We are interested in possible collective excitations, 
the role of shell structure, and in general in the factors de¬ 
termining the amount of strength. As in previous work, 
we calculate the properties of dipole excitations using the 
QRPA and the Gogny D1S interaction. Furthermore, we 
use a continuum-RPA model (CRPA) and the Skyrme 
functional SLy4, to examine the threshold strength. We 


will compare our results with available data and other 
theoretical results. 

Our dipole-mode nomenclature should be clarified at 
this point. In the following we refer repeatedly to two 
types of collective dipole excitations. One is an ordi¬ 
nary surface dipole vibration , which corresponds to the 
isoscalar low-energy dipole (IS-LED) mode studied in 
Refs. [15-19]. We will therefore denote it as IS-LED. 
The other is driven by neutron excess and we will refer 
to it as a neutron mode. The term neutron-skin mode or 
oscillation we reserve explicitly for the oscillation of the 
neutron skin as a spatially defined entity. Some authors 
call that the pygmy resonance. In this work, however, we 
use the term pygmy to refer to low-energy El transitions 
of any physical origin, as long as they are not considered 
fragments of the giant dipole resonance. 

The plan of this paper is as follows. In Sec. II we intro¬ 
duce briefly the theoretical models we use and define the 
quantities of interest. In Sec. Ill we present and analyse 
our results. In particular, in Sec. Ill A we examine the 
isovector and isoscalar response of Ni isotopes and iden¬ 
tify interesting low-energy structures to be studied next. 
In Sec. Ill B we examine the evolution of those structures 
with respect to the neutron number. Having identified 
two potentially interesting modes, in Sec. Ill C we an¬ 
alyze them through their transition densities and wave- 
functions. In Sec. Ill D we calculate the total low-energy 
strength and the dipole polarizability and compare our 
results with existing data. A critical discussion of our 
findings is presented in Sec. IV. We focus on the rele¬ 
vance of shell structure and the phenomenon of isospin 
splitting and we comment on how our results fit into the 
ongoing discourse about the relation of dipole states with 
the softness of the symmetry energy. We summarize our 
conclusions in Sec. V. 


II. THEORY BACKGROUND 
A. Theoretical models 

In this work we employ two models based on the 
random-phase approximation. The models have been dis¬ 
cussed in detail elsewhere. Therefore, we will now present 
them only briefly and refer the reader to respective pub¬ 
lications for more details. 

We use mainly the self-consistent quasi-particle 
random-phase approximation (QRPA) with the Gogny 
D1S interaction [32], which has yielded already interest¬ 
ing results in our previous studies of dipole strength [16, 
17, 19]. The QRPA is formulated in the canonical basis 
of the Hartree-Fock-Bogolyubov model (HFB) describ¬ 
ing the ground state. The same Hamiltonian is used in 
the description of the ground states and in the ph and 
pp channels in QRPA, leading to an excellent degree of 
self-consistency. The theory and implementation are pre¬ 
sented in detail in Ref. [33]. The reader may also consult 
Ref. [17] for a concise presentation. A harmonic oscillator 
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basis of 15 shells is used in the present work. The oscilla¬ 
tor length parameter is chosen for each isotope such that 
the energy of the HFB ground state is minimized. 

Our QRPA model in the harmonic-oscillator basis lacks 
a proper treatment of the continuum. For this reason 
we also explore the continuum random-phase approxima¬ 
tion (CCRPA) for Skyrme interactions. We employ the 
SLy4 parameterization [34], because it predicts similar 
nuclear-matter properties of relevance as the D1S model, 
namely the symmetry energy and the nucleon effective 
mass. Both models are asy-soft. They differ mainly in the 
Thomas-Reiclie-Kuhn (TRK) enhancement factor, which 
influences the energy of the GDR, but not strongly the 
pygmy-strength properties or the neutron skin [5, 11, 35]. 

The CRPA formalism makes use of Green’s functions. 
The ground state is obtained with the Skyrme-Hartree- 
Fock model [36]. The unperturbed ph propagators and 
the residual interaction are evaluated in coordinate space. 
Convergence issues with respect to the model space are 
thus automatically absent. Appropriate boundary con¬ 
ditions account for the possibility of particle emission 
when the particle energy is positive. The escape effect 
can be turned off by imposing box boundary conditions, 
namely that all wave functions vanish at some distance 
from the nucleus. The model is the same as that used in 
Ref. [37] and the reader is referred there for further infor¬ 
mation. The CRPA does not include pairing and we will 
use it only for a few isotopes which afford a description 
as closed-shell configurations. We typically use a radial 
mesh of 18 frn and a step size of 0.08 frn. 


B. Quantities of interest 


The transition strength distribution corresponding to 
an external field Q is given, in QRPA, by 

S{E) = Y j B{E v ) 5{E-E v ), (1) 

v> 0 


where the transition strength of the v— th eigenvalue 
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is given in terms of the 2-quasiparticle X, Y amplitudes 
and the occupation probability factors u,v, 

bqq' = i X qq' ~ Y qq')(VqUq' ~ V q > U q ) (<?| Q\q') . (3) 

The minus sign follows from J 7r = 1~. In CRPA the 
strength distribution is determined by the polarization 
diagram via the Green’s function 

S{E) = -~qY, f <Prd\'f t (r)G tv {ry-,E)f tl {r'), 

( 4 ) 


where ft.{r) denotes the spatial and isospin dependence 
of the external field. 

The isoscalar (IS) and isovector (IV) electric dipole 
response is determined in the two models, in the long- 
wavelength limit, by the following respective operators, 
which are corrected for spurious center-of-mass effects: 

Qis = f r i - t(r 2 )r i ') V^hWG), (5) 

i=1 ' ' 

and 

N . Z ^, 

Qei = e—^2r p Y 1M (rp) - e—^2r n Y 1M (r n ). ( 6 ) 

p— 1 n— 1 

As far as intrinsic excitations are concerned, the operator 
Qei is equivalent to the uncorrected El operator and to 
the isovector dipole operator 

z A 

Qei o e^V p Yi M (Ap) <J=> § X] r 3* )r ViM(Aj), 

P= 1 2=1 

in an obvious notation. When we speak of El strength 
we will refer to the electric-dipole excitation strength 
B{E1) f corresponding to the operator Qei and we shall 
use IS strength when we refer to the operator Qis- 

The QRPA transition strengths obtained with the cor¬ 
rected operators (5) and (6) are the same as with the un¬ 
corrected forms, except of course for the spurious state. 
The spurious state appears always very close to zero en¬ 
ergy. 

Our CRPA implementation is not fully self-consistent, 
in that spin-dependent operators and the Coulomb term 
are omitted from the residual interaction. Therefore, we 
always scale the residual interaction by a factor such that 
the spurious state appears close to zero energy. The scal¬ 
ing factor ranges from 1.06 for 68 Ni, to 1.11 for 48 Ni, to 
the large values of 1.15 for 78 Ni and 1.22 for 56 Ni. The 
use of corrected operators is imperative. 

The discrete energy spectra obtained with QRPA can 
be smoothed with a Lorenzian of a given width, for pre¬ 
sentation purposes. In CRPA the spectra are continu¬ 
ous and include the escape width. In order to visual¬ 
ize very narrow or bound structures, an artificial width 
can be introduced by adding an imaginary constant to 
the particle-hole energy. In our presentation of results, 
when we say that our QRPA or CRPA results have been 
smoothed by means of a given width parameter, we will 
imply the width at half maximum. 

The dipole polarizability 

he I (Tabs , , > 

= (7) 

which is proportional to the inverse-energy-weighted sum 
of El strength, has been found to correlate with the 
neutron-skin thickness and the density-dependence of the 
symmetry energy [5]. It has been measured in 208 Pb [7], 
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120 Sn [8], and in 68 Ni [38]. We will calculate this quan¬ 
tity along the Ni chain. In order to eliminate a trivial 
dependence on the mass number, we define a reference, 
collective value a'g' 11 , as the polarizability corresponding 
to the giant dipole resonance if it was a single peak with 
the empirical energy [39] 

E gdr = (31.2A” 1/3 + 20.6A" 1 / 6 ) MeV 

and exhausting 100% of the classical TRK sum rule. This 
gives us 1 

coll he 60 N ZMeVrnb 
“ D = ^ (8) 

neglecting the TRK enhancement factor. 

The transition density, between the ground state and 
an excited state, determines how the system couples to 
the external field at the given excitation energy and pro¬ 
vides insight into the microscopic structure of the excited 
state. In QRPA, solved in a spherical basis, the radial 
part of the transition density is given by 

Spt = ( X qq' - Y qq')( V q U q' ~ V q' Uq)R q (r)R q > (r ), 

q<q',tq-t'q 

where R q (r) is the radial wavefunction of the canonical 
state and t q refers to its isospin (proton or neutron). 

In our CRPA model, transition densities are evaluated 
as follows. In the vicinity of a resonance, at energy E, the 
Green’s function is proportional to a product of transition 
densities [41]. In particular, 

$sG t , t '(r,r'-,E) oc 5p t {r-,E)5p t >{r';E), (9) 

where we have used spherical symmetry to separate out 
the angular dependencies, and t denotes a proton or a 
neutron. Folding the above quantity with the external 
field, we get a function 

Sp { t f \r-,E) = E / %G ttl {ry-,E)f t ,{r')dr', (10) 

t' 

which is proportional to 8pt{r\E), if Eq. (9) holds. We 
observe that if we fold this with the external field / t (r), 
we obtain 7r times the transition strength (per energy 
unit). The above factorization is not always perfect, but 
in any case the above procedure helps us obtain a quan¬ 
tity, loosely called a transition density, which determines 
how the nuclear density couples to the external field at a 
given energy. It may depend on / t (r), therefore, we use 
the subscript /. 


1 After submitting this manuscript we learnt that a similar refer¬ 
ence value has been defined by S.Ebata et al. for related stud¬ 
ies [40]. 


C. Assessing collectivity 

Different complementary ways can be used to assess 
the collectivity of a QRPA eigenstate. One of them is 
the amount of 2-quasiparticle configurations which con¬ 
tribute considerably to the norm or to the transition 
strength (summation in Eq. (2)) of a given QRPA eigen¬ 
value. In the case of dipole states, care must be exer¬ 
cised when using the latter criterion, as the results de¬ 
pend on the excitation operator used (i.e., whether it is 
corrected or not) [19]. Therefore, in addition, we will 
employ single-particle units, in order to assess the col¬ 
lectivity of a given state. For El transitions the usual 
single-particle unit must be corrected for the center-of- 
mass motion, by multiplying the usual estimate [42] by 
the effective-charge factor squared, namely ( N/A ) 2 for 
proton states and ( Z/A ) 2 for neutron states [39]. We 
choose to use the single-neutron unit 

B s . n .(El) = 0.06445A 2/3 (Z/A) 2 e 2 fm 2 . (11) 


III. RESULTS, ANALYSIS, AND COMPARISON 
WITH DATA 

A. Main features of response 

First we compare the general features of the response, 
in particular giant resonances, between the two models 
and with some data. The only isotopes which can be de¬ 
scribed using closed subshells and, therefore, lend them¬ 
selves to calculations with both models are those with 
N = 20,28,40,50, all of which are unstable. Relevant 
data are available for stable isotopes and for 68 Ni. We 
will, therefore, use 58 ’ 68 Ni for a first comparison of our 
results with data. 

The nucleon separation energies for the N = 
20,28,40,50 isotopes are calculated within the D1S 
model to be 0.41,7.31,8.59,5.59 MeV, respectively. For 
N = 28,40 the values given in the AME2012 evalua¬ 
tion [43] are 7.167, 7.793 MeV, respectively. 

In Fig. 1 we show the calculated photoabsorption cross 
sections for the above isotopes. A large smoothing fac¬ 
tor has been used. The giant dipole resonance (GDR) 
is visible in all cases, with some degree of fragmenta¬ 
tion. The GDR energies predicted by the two models 
are visibly different, by a few MeV. The difference can 
be attributed to the different TRK enhancement factor 
(0.6 for D1S, 0.25 for SLy4). An inspection of GDR data 
indicates that the QRPA+D1S model reproduces more 
correctly the higher-energy part and overall the mean 
energy of the GDR, while the CRPA+SLy4 model pre¬ 
dicts rather correctly the main first peak of the GDR. 
We infer this from Fig. 1(d) and from a juxtaposition of 
Figs. l(b,c) in the expectation that the GDR properties 
of 56 Ni and 58 Ni must be quite similar. Both models can 
be considered acceptable in this sense. Finally, they both 
appear to underestimate the fragmentation of the reso- 
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nance. Indeed, the fragmentation and spreading of the 
GDR can be enhanced by effects beyond (Q)RPA, like 
collisional damping [44], especially in heavy nuclei, and 
clustering [45], especially in light nuclei. In medium-mass 
nuclei both mechanisms may contribute. 

In Fig. 2 we show the calculated isoscalar strength dis¬ 
tributions. In 58 Ni IS dipole strength has been detected 
both below [49] and above [50] particle emission thresh¬ 
old and the major structures are indicated in Fig. 2(c) 
with arrows. 

In all 5 isotopes, the broad structure at high energy, 
around 30 MeV, corresponds to the IS GDR, or dipole 
compression mode. We have confirmed that the pro¬ 
ton and neutron transition densities show the typical 
behavior for this mode [39], namely they are in phase 
and have a node. The smooth distribution given by the 
CRPA-|-SLy4 model in this energy region is not the result 
of a large smoothing factor, but of the genuinely large es¬ 
cape width, reproduced through the appropriate bound¬ 
ary conditions. The use of box boundary conditions (not 
shown) produces a structured distribution similar to the 
QRPA+D1S case. In 58 Ni, the high-lying peak was ob¬ 
served at 34 MeV, while the calculated response func¬ 
tion peaks around 30 MeV. In heavier nuclei theoretical 
models tend to predict a higher energy for the IS GDR 
than extracted experimentally [51]. The reason for these 
discrepancies is likely a combination of factors, not re¬ 
stricted to the compression modulus, since discrepancies 
remain when the monopole compression mode is correctly 
described. 

At least two more structures, on which we shall remark 
next, are visible at lower energies: around 10 MeV and 
around 15 — 20 MeV, depending on the isotope and the 
interaction. 

The lowest IS peak, typically close to 10 MeV, is the 
IS-LED mode (or surface mode), and we analyze it in 
Sec. Ill C. In 58 Ni IS dipole strength has been detected 
in this energy region, as indicated in Fig. 2(c). In 48 Ni, 
this peak would lie within a broad continuum of proton- 
dominated transitions, as we infer from the CRPA+SLy4 
results of Fig. 2(a) and from the similar results analyzed 
in [37]. The large width obtained with the CRPA+SLy4 
model for the low-energy strength distribution of this nu¬ 
cleus reflects the proton escape width. 

The structure at 15 — 20 MeV is no less interesting. 
Such an IS dipole peak, just below the IV GDR, has 
been detected in a variety of nuclei [51, 52], including 
in 58 Ni at 17.4 MeV. Its physical interpretation remains 
unclear [53-55]. A detailed study of this structure lies 
beyond the scope of the present work. In the following we 
will focus on the structures closer to and below 10 MeV, 
relevant to the physics of pygmy resonances. 

We may already observe that both models predict in¬ 
creased El strength in 68 Ni and ,8 Ni around 10 MeV-cf. 
Fig. 1(d),(e). Next we will identify transitions of poten¬ 
tial interest which give rise to this and other effects. 


B. Systematics with neutron number 

Fig. 3 gives an overview of the response of all isotopes 
calculated with the QRPA+D1S method. The IS and El 
transition strengths are indicated by green disks and open 
black circles, respectively. The area of a disk or circle is 
propotional to the strength. The basic features of the 
spectrum and how they evolve with neutron number are 
clearly visible. The large circles in the energy region of 
20 ± 3 MeV correspond to the GDR, which appears frag¬ 
mented into two or more structures. The lower-energy 
part of the GDR becomes strongly isoscalar in the more 
asymmetric isotopes, thus showing a mixed isospin char¬ 
acter. The disks at high energy, above 25 MeV, corre¬ 
spond to the high-lying compression mode. 

The low-energy spectrum in the IS channel is domi¬ 
nated, in all isotopes, by at least one strong transition 
close to 10 MeV. We may interpret the lowest-energy 
peak, in all isotopes with N = 20 — 50, as a mode of sim¬ 
ilar nature as the IS-LED studied in Refs. [16-19]. We 
study its properties in Sec. IIIC. Even in 80,82 Ni, the 
states visible at 9.5 and 9.7 MeV are candidates for such 
a mode. 

As neutrons are added to the N = Z nucleus 56 Ni, this 
transition appears to aquire more and more El strength. 
When we reach 62 Ni we observe a splitting in the low- 
energy IS strength into two promiment transitions. Both 
appear to be of mixed isospin character, with the lower 
one remaining more strongly IS. The second structure 
presents an interesting candidate for an exotic mode gen¬ 
erated by neutron excess. We will also analyse it in 
Sec. IIIC. 

In the proton-rich isotopes the IS-LED carries a very 
small amount of El strength. Based on our results with 
CRPA+SLy4 for 48 Ni and those of Ref. [37], we expect it 
to lie in the continuum and mix with proton transitions 
in most of these isotopes. Finally, beyond the shell clo¬ 
sure of N = 50 a large amount of IS and El strength is 
obtained at lower energies and it all is expected to lie in 
the continuum (see also the results in Refs. [56, 57]) along 
with the IS-LED we identify in 80,82 Ni. The nucleon sepa¬ 
ration energies in 80 - 82 . 84 . 86 ]\[i ; calculated within the D1S 
model, equal 2.05,1.88,1.84, 0.71 MeV, respectively. 

C. How many modes? 

We now examine the two types of low-energy states 
which we have identified as an isoscalar mode (IS-LED) 
excited in all isotopes, and a neutron mode in 62 ~' 6 Ni. 

The coherent IS-LED mode of 56 Ni has been analyzed 
in Ref. [16], along with that of other N = Z nuclei. The 
corresponding proton and neutron transition densities are 
very similar to each other and have a node, fitting well in 
the picture of the surface dipole mode of ordinary nuclei 
discussed in Ref. [15]. Let us now see to what extent such 
a picture persists in the other Ni isotopes. 

In Fig. 4 we show the proton and neutron transition 
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FIG. 1. (Color online) Calculated photoabsorption cross section, for the indicated interactions and isotopes. The calculated 
cross sections are smoothed with a Lorenzian of width equal to 2 MeV. For 58 Ni (c) the dashed curve corresponds to a smoothing 
width of 5 MeV. The data plotted in (c) and (d) for 58 Ni and 68 Ni are from Refs. [46] and [38, 47], respectively. The arrows in 
(d) indicate peaks identified in Refs. [38] (black) and [48] (gray). 
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FIG. 2. (Color online) Calculated isoscalar strength distribution, for the indicated interactions and isotopes. The distributions 
are smoothed with a Lorenzian of width equal to 2 MeV. For 58 Ni (c) the dashed curve corresponds to a smoothing width of 
5 MeV. The arrows show peaks detected below 10 MeV in (a, a'jo) experiments [49] and mean energies of gaussian fits to the 
IS strength distribution extracted from alpha scattering at small angles [50]. 


densities for the IS-LED in all even isotopes with A = 
48 — 82, calculated within the QRPA+D1S model. For 
A = 48 — 78 that is the lowest-lying IS state in Fig. 3. For 
80 ’ 82 Ni it is the IS state at 9.5 and 9.7 MeV, respectively. 
We notice that in all isotopes both transition densities 
are characterized by a node and they are in phase. 

The transition densities change very smoothly from 
48 Ni to 56 Ni and from 70 Ni to 82 Ni, and quite smoothly 
in between the two domains. In the neutron-deficient do¬ 
main and including 56 Ni the similarities between the pro¬ 
ton and neutron transition densities and in all isotopes 
are obvious. All these states can, therefore, be confirmed 
as of the usual IS-LED character. 

For A > 60 the neutron transition densities are 
stronger than the proton counterparts and somewhat dis¬ 
similar to them. Nevertheless, the surface protons con¬ 
tribute even in the very neutron-ricli isotopes. The jux¬ 
taposition with the plotted lowest-lying state in 80 Ni un¬ 
descores this point. Therefore, to some extent we may 
also consider them as IS-LED states. 

The observed discontinuities, with respect to the neu¬ 
tron number, especially visible in the neutron transition 
densities, can be attributed to shell structure, namely the 
filling of the neutron pf shell and finally the 0 <? g / 2 state. 
In reality, some of the discontinuities may be smoothened 
out by correlations beyond HFB and QRPA. The occu¬ 


pation probabilities of the above key orbitals, calculated 
within the HFB method, are shown in Fig. 5. From 60 Ni 
to 68 Ni the almost-degenerate (in this model) lpi/ 2 , 0 / 5 / 2 
states fill up at a similar pace, while the lpi / 2 state is al¬ 
most full and the 0g 9 / 2 empty. Beyond 68 Ni the strongest 
difference is the filling up of the 0<? g / 2 orbital. 

In Fig. 6 we show the transition densities of the 
neutron-mode candidates in 64 ’ 68 , 72 ’ 76 Ni. For compar¬ 
ison, the IS-LED transition densities are also plotted 
with thinner lines. There are similarities between the 
two modes, in particular the presence of nodes and the 
non-negligible contribution of surface protons in most 
cases. We do notice, however, that the proton contri¬ 
bution to the neutron-mode candidate diminishes in the 
more neutron-rich isotopes - see Fig. 6 (c),(d). Further¬ 
more, this mode is predicted to decouple from the IS-LED 
and to develop only in neutron-rich Ni isotopes. For this 
reason we identify this as a separate mode, driven by 
neutron excess and assisted by shell effects. In particu¬ 
lar, 62 Ni is the lightest isotope for which the lp \/2 and 
O/ 5/2 states are significantly occupied, activating new 2 qp 
configurations. An analysis of the wavefunctions shows 
that the 0 / 7/ 2 —> 0g 9 /2 and the lpOf —► 2sld transi¬ 
tions are coupled into the IS-LED for the lighter isotopes, 
while for 62 Ni and beyond they split from each other and 
lpOf 2sld attract and couple to the newly available 
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FIG. 3. (Color online) For the given isotopes and excitation 
energies, the IS and El transition strengths are indicated by 
green disks and black circles, respectively. The area of a disk 
or circle is propotional to the strength. 


excitations of 0 gg /2 neutrons at somewhat higher ener¬ 
gies. 

The lowest-lying states in 80-86 Ni (below 8 MeV in 
Fig. 3) can be attributed to loosely-bound neutrons in 
spatially extended orbitals of the sdg shell. These con¬ 
stitute yet another type of transition, similar to the soft 
modes of halo nuclei. The neutron transition densities 
extend well beyond 8 fm, while the proton transition 
densities remain localised and peak between 4 and 5 fm 
(surface-peaked). This is exemplified by the transition 
density of the 7.3 MeV state of 80 Ni, shown in Fig. 6. 
The strongly isoscalar and very weakly isovector transi¬ 
tions just above 10 MeV resemble the IS-LED, in that the 
proton transition densities have a node. We may there¬ 
fore distinguish two groups of low-lying strongly isoscalar 
states beyond 80 Ni, depending on how strongly isovector 
they are. Such broad structures might be observable, 
though not as individual transitions, because they are 
unbound. 

We now examine whether the IS-LED and the neutron 
mode in the lighter isotopes (A < 78) can be considered 
collective. The IS-LED is found strong in the IS chan¬ 
nel, exhausting a few percentage points of the energy- 
weighted sum, and it energetically lies below the several 
lHui configurations contributing to its IS strength. We 


A = 48 (thickest) - 56 (thinnest) 
A = 58 (thickest) - 62 (thinnest) 
A = 64 (thickest) - 68 (thinnest) 
A = 70 (thickest) - 78 (thinnest) 
A=80, 7.3 MeV 
A=80, 9.5 MeV 



FIG. 4. (Color online) (a) Proton and (b) neutron transition 
densities of the IS-LED mode in even Ni isotopes, calculated 
with the QRPA+D1S model. For comparison, the transition 
densities of the lowest dipole mode in 80 Ni are also shown. 
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FIG. 5. (Color online) Calculated occupation probabilities of 
the neutron 0flp0gg/2 states in 56 ~ 78 Ni. 
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may, therefore, consider it as coherent-see also Ref. [16]. 
The neutron-mode candidate is quite strong in the IS 
channel around 68 Ni and lies below most IHuj configura¬ 
tions. Its El strength exceeds one single-neutron value, 
as we will see in Sec. HID when we discuss Fig. 9. The 
above are signs of collectivity. 
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FIG. 6. (Color online) Proton and neutron transition densities of the neutron-mode (or pygmy) candidate in 64,68,,2 ’ 76 Ni, 
calculated with the QRPA+D1S model. For comparison, the transition densities of the IS-LED mode (from Fig. 4) are also 
shown, plotted with thinner and paler lines. 
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FIG. 7. (Color online) The contributions of 2-quasiparticle 
configurations to the transition matrix element of the IS-LED 
and the 10.6MeV mode in 68 Ni, calculated in the QRPA+D1S 
model, (a) for the IV operator, Eq. (6), and (b) for the IS 
operator, Eq. (5). 
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It is worth examining the wavefunctions of this type 
of mode. We will look at the matrix element terms b u qq ,, 
Eq. (3), which determine the contribution of each config¬ 
uration to the strength of the 10.6—MeV mode in 68 Ni. 
The above terms for the IV and IS operators are plotted 
in Fig. 7(a) and (b), respectively, as a function of the 
2-quasiparticle energy. For comparison, the same quan¬ 


tities are shown for the IS-LED. We use the corrected 
operators, Eqs. (5), (6). Notice that if we considered 
the uncorrected IV operator in Fig. 7(a), only the pro¬ 
ton contributions would appear in Fig. 7(a), multiplied 
by a factor A/Z, though the final result for the strength 
would be the same. Similarly, use of the uncorrected IS 
operator would alter the results in Fig. 7(b). We stress, 
therefore, that there is some ambiguity in analyses such 
as the following. 

The prominent proton and neutron contributions to 
the El strength at 12.8 and 12.7 MeV, respectively, in 
both modes come from the f ^/2 —> < 79/2 configurations, 
which are not spatially extended. There are consider¬ 
able destructive contributions from many other config¬ 
urations. The neutron configurations at 12.2 and 13.9, 
namely IP 3/2 -A lds /2 and lsi /2 -A ICI 3 / 2 , are the ones 
which contribute more strongly to the neutron mode, 
than the IS-LED. They are low —l configurations and 
therefore spatially extended. On the other hand, the 
v2s\/2 —> lpi /2 configuration at 12.66 MeV contributes 
only to the IS-LED. 

In Fig. 7(b) we see that the IS contributions of the 
various lhu) configurations are mostly coherent in both 
modes. We conclude from the above that both can be 
considered coherent IS modes, made up of fragments of 
many ltuo configurations. In this respect our results 
are in agreement with those of Ref. [10], obtained using 
Skyrme forces, at least as regards the sole mode analyzed 
there. 

Based on the above analysis, and on the much higher 
El strength of the second mode at 10.6 MeV, we may 
consider the latter as of a different kind, namely an exci¬ 
tation of neutrons on the nuclear surface. The presence of 
the vfr /2 -A ffg /2 configuration appears essential in gen¬ 
erating coherence for this mode. Therefore the strength 
diminishes in 78 Ni, where the neutron 0<?g/2 state is filled. 

In Fig. 8 we show the IS and El strength distribu¬ 
tion of 68 Ni and 78 Ni in this energy region, calculated 
with both models. We notice that the bimodal struc¬ 
ture of the response of 68 Ni is present also in our results 
with the CRPA+SLy4 model. Noteworthy are the sig- 
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(Color online) The IS and El strength distribu¬ 
tions, for 68 Ni and 78 Ni, calculated with the QRPA+DIS and 
CRPA+SLy4. The response functions are smoothed with a 
width of 0.4 MeV. The data in (c) are from Ref. [38, 47). 


nificant contributions from the non-resonant continuum 
in both channels in the region up to around 11 MeV. 
The effect is predicted to be especially important in ' s Ni. 
Such continuum transitions, first studied in the 1990’s 
(see, e.g., [56, 57]), are likely critical for describing cor¬ 
rectly the supratlireshold strength beyond 68 Ni. In fact, 
in the case of 68 Ni, if both resonances, namely the IS- 
LED and the neutron mode, are bound transitions, the 
QRPA+DIS cannot explain the observed strength. The 
same can be said about other configuration-space calcu¬ 
lations, employing the SLy4 or other Skyrrne function¬ 
als [58]. In general, the energies of the unbound single¬ 
particle states, and therefore the non-collective single¬ 
particle spectrum in the continuum, converge very slowly 
with the harmonic-oscillator basis size. In actual calcula¬ 
tions they can never be considered converged. If the par¬ 
ticle threshold is much lower than the GDR, this short¬ 
coming becomes relevant for pygmy strength. 

It is noteworthy that the CRPA+SLy4 model predicts 
increased El strength also around 12 — 13 MeV, com¬ 
pared to the QRPA+DIS model, in the form of additional 
peaks. Peaks in this energy region are obtained for all 
four isotopes studied with the CRPA-SLy4 model, as is 
evident in Fig. 1. As the most possible reason for the ad¬ 
ditional CRPA strength in that energy region we find the 
overall energetic displacement of one GDR strength dis¬ 
tribution with respect to the other. The proton and neu¬ 
tron transition densities r 2 5p l ^ n (r) of the CRPA peaks 
in question, for N > Z , typically have prominent peaks 
at a value of r somewhat larger than the mean-square 
radius and they are out of phase at that r—value. 2 Then 



FIG. 9. (Color online) Sum of El strength for the indicated 
energy cuttoffs as a function of mass number. For 58,60 Ni, 
data shown are from Refs. [59, 60]. The data points for A = 68 
are: an estimate assuming a single peak with the reported 
mean energy and energy-weighted strength [38] (lower point) 
and the actual sum of strength measured between 7.8 and 
10 MeV [47] (upper point). An earlier measurement of the 
energy-weighted summed strength [48, 61] gives the estimate 
l.l±0.3e 2 fm 2 . The single-neutron unit, B s n (E 1), is given by 
Eq. (11). 


the peaks could be considered as fragments of the GDR. 
They are, however, embedded in a continuum of transi¬ 
tions dominated by surface neutrons. We note that the 
existing data for 5S Ni (Fig. 1(b)) and 60 Ni (not shown) 
do not present a prominent structure on the lower-energy 
tail of the GDR. We obtained similar results for the GDR 
and for these peaks using the SkM* functional instead of 
the SLy4. A similar fragmentation pattern is observed 
even with the QRPA+DIS model at somewhat higher 
energies. 


D. El transition strength and polarizability 

The El strength of the IS-LED and the summed El 
strength for different energetic cutoffs as a function of 
mass number is shown in Fig. 9. First we discuss our 
results with the QRPA+DIS model. We observe that 
the El strength of the IS-LED rises up to A = 60 and 
then drops again, as a result of the splitting already seen 
above. An interesting result is that the summed strength 
up to 12 MeV or similar value does not rise monotoni- 
cally, but it peaks close to 68 Ni. Similarly, it has been 
observed that the low-energy El strength does not rise 


2 In 48 Ni (and to a large extent in 78 Ni) the surface protons (neu¬ 


trons) play a dominant role. 
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monotonically between 40 Ca and 48 Ca [62]. Clearly, the 
neutron excess is not the only mechanism determining 
the El strength. In the context of our calculations, the 
explanation for the above behavior must lie in the shell 
structure. As we already mentioned, between 68 Ni and 
,s Ni the 0<?9/2 state fills up, blocking transitions from 
the O/ 7/2 state. Although transitions from the 0gg/2 are 
available for N > 40, we found that they contribute little 
strength in this region. The above anticorrelation be¬ 
tween the strength and the neutron number in neutron- 
rich Ni isotopes was found also in Ref. [63] and attributed 
to the filling of the 0 g 9 / 2 orbital. 

The strength of the neutron mode in 68 Ni, analyzed 
in Sec. IIIC, exceeds the single-particle estimate. The 
summed El strength in this region was found large in 
Sn isotopes too [19], but in those cases it is generated by 
individual weak transitions. 

Given our past results and how they compare with 
data [17-19], we expect the summed strength up to 
12 MeV to correspond to a summed experimental 
strength up to, roughly, 9 MeV. For the stable isotopes 
this would mean an excellent agreement with the data 
shown in Fig. 9. In 68 Ni, pygmy strength was found in 
the continuum up to 12 MeV. Notice that, again based 
on the known trends of the model, we expect the first 
IS peak to be bound, i.e., below the threshold energy of 
7.8 MeV. (For 72 Ni and heavier isotopes, for which the 
neutron threshold is below 7 MeV, practically all strength 
is likely in the continuum.) It follows that our results are 
marginally compatible with the data in the case of 68 Ni. 

Now we turn to the El strength obtained with the 
CRPA+SLy4 model, also shown in Fig. 9, for various 
summation cutoffs. We notice that the predictions of 
this model agree with the data for 68 Ni better than those 
of the QRPA+D1S model. A much more dramatic differ¬ 
ence between the two models is seen in 78 Ni (cf. Fig. 8). 
Threshold transitions described efficiently by the CRPA 
model likely become increasingly important beyond 68 Ni. 

Finally we present results for the dipole polarizability 
an of the Ni isotopes. In Fig. 10 we plot an relative to the 
reference value a^ 1 defined in Eq. (8), to eliminate trivial 
dependencies on the mass number. The polarizability as 
a function of neutron number varies in a smooth manner 
between shell closures, at which points, namely 56 Ni and 
78 Ni, it presents kinks. Obviously this quantity is influ¬ 
enced by shell structure. Beyond 78 Ni it rises rapidly. 
This behavior can be attributed to the strong low-lying 
transitions observed in Fig. 3. Indeed, if we calculate an 
without taking into account the lower-lying states, the 
function cld(A) / a c £ ll (A) becomes smoother and follows 
closely the linear trend, as evident in Fig. 10. 

We notice that the two models overestimate the polar¬ 
izability of 68 Ni, even though the functionals are asy-soft 
and underestimate the strength of the observed pygmy 
state. We elaborate on this contradiction at the end of 
Sec. IV B. 

In Fig. 10 the neutron or proton-skin thickness, defined 
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FIG. 10. (Color online) The electric dipole polarizability ac 
of Ni isotopes (left axis) calculated with both models as in¬ 
dicated and the difference of the proton and neutron root- 
mean-square radii (right axis) calculated with the HFB+D1S 
model. The polarizability is given relative to the estimate o^ 11 
for a single collective mode. The radii difference is divided by 
A 173 . The measurement of a_o in 68 Ni is from Ref. [38]. 


via the root-mean-square radii 


I E tj Rp 




is also shown as a function of A. In particular, we plot the 
quantity \R n — i? p |/A 173 , which represents a skin thick¬ 
ness relative to the nuclear radius. It is calculated in 
the HFB model with the Gogny D1S interaction. The 
value obtained for the neutron-skin thickness in 68 Ni is 
0.15 frn. The proton and neutron radii are found al¬ 
most equal in 58 Ni, while for A < 56, i.e., Z < N, 
a proton skin develops. We notice that the (proton or 
neutron) skin thickness increases monotonically with the 
proton-neutron asymmetry, except around N = Z. The 
effect of neutron deficiency is stronger when the contin¬ 
uum is propertly treated (CRPA+SLy4). Importantly, 
the fact that the relative polarizability decreases beyond 
'°Ni shows that it does not correlate with the relative 
neutron-skin thickness, but it is strongly influenced by 
shell effects. A possible correlation emerges beyond the 
N = 82 shell closure - see also Ref. [64]. We conclude 
from the above that the relation between the proton- 
neutron asymmetry and the polarizability, along an iso¬ 
topic chain, is not completely trivial. 


IV. CRITICAL DISCUSSION IN VIEW OF 
OTHER EXISTING CALCULATIONS 

A. Ca, Sn, isospin splitting, and the importance of 
shell effects 

We now compare results for the three isotopic chains, 
Ca [17, 18], Ni (this work), and Sn [19]. 
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We begin with the following general observation, which 
will help us interpret the isospin structure of pygmy 
strength and especially the so-called isospin splitting in 
heavy stable nuclei. By the latter term we mean the 
phenomenon whereby El transitions are observed over a 
range of excitation energies below particle threshold, of 
which only the lowest-lying respond to an IS field [1]. 

First we observe, in the compilation of Ref. [15], that 
the energy of the isoscalar segment, typically 6—7 MeV, 
varies weakly with mass number in comparison, e.g., 
to the shell energy 1 hio = 41A _1 / 3 MeV. Pygmy dipole 
strength, according to non-relativistic models, is ex¬ 
pected at around 1 hut [12] and can be attributed to frag¬ 
ments of 1 hui transitions, rather than a collective mech¬ 
anism [3, 9-11]. In medium-mass stable nuclei, such as 
stable Ca and Ni isotopes, the IS and pygmy (1 Hco) tran¬ 
sitions generally lie far apart, because of the above A 
dependence. In heavier stable nuclei, at least the lower- 
energy tail of the pygmy strength distribution can ap¬ 
proach and overlap with the IS state. 

Indeed, a noteworthy difference observed in the exper¬ 
imental data on stable nuclei is that just a few weak 
IV transitions are seen on either side of the IS-LED in 
48 Ca [18], while the IS-LED candidate lies on the lower 
end of the, less weak, El spectrum in stable Sn iso¬ 
topes [17, 65] and other heavy nuclei [1]. A plausible 
explanation for the so-called isospin splitting of pygmy 
strength in heavy stable nuclei is then the energetic prox¬ 
imity of two resonance-like structures: the IS-LED and 
the tail of the higher-lying pygmy mode (whether collec¬ 
tive or not). 

We conclude from the above that the isospin character 
of the dipole strength distribution at low energies de¬ 
pends on the mass region. 

Regarding the isospin structure of the El states in sta¬ 
ble Ni isotopes, below threshold, we expect that the IS 
and (weak) IV strength distributions will mostly coexist 
without a clear energetic splitting. Our expectation is 
based on the extended IS strength distribution detected 
in 58 Ni [49]. The situation would resemble that in 48 Ca, 
but with a fragmented IS strength distribution. We can¬ 
not conclude on the matter based on our calculations, 
since QRPA predicts the IS-LED always to be the lowest- 
lying dipole mode (see Ref. [18] regarding 48 Ca) and does 
not include fragmentation mechanisms like phonon cou¬ 
pling [21-24], 

Next we turn to a different observation. In the Ca and 
Sn isotopic chains we found that the El strength of the 
IS-LED state does not increase monotonically as we add 
only a few neutrons to the N = Z isotope, but it reaches 
a minimum for N slightly larger than Z. Since the El 
strength of this IS mode is a result of local dissimilar¬ 
ities between the proton and neutron transition densi¬ 
ties [16], reflecting, in turn, dissimilarities in the ground- 
state densities, it appears easy to explain its behavior as 
an effect of the Coulomb interaction: For N = Z the pro¬ 
ton density distribution should be slightly more extended 
than the neutron density distribution, leading to this lo¬ 


cal dissimilarity. Addition of a few neutrons reverses the 
situation, resulting in minimal El transition strength in 
between, at some small value of N — Z. 

This interpretation, which relies on a senri-classical im¬ 
age of density profiles, breaks down in the Ni chain, where 
minimal El strength is obtained for N slightly smaller 
than Z , signifying a competition between Coulomb and 
asymmetry effects, on the one hand, and shell effects on 
the other. Indeed, the Ni shell structure tends to pro¬ 
duce the opposite effect to asymmetry: Addition of neu¬ 
trons to 56 Ni means filling up lower-1' states, which are 
extended to the nuclear surface, immediately initiating 
what resembles a neutron-skin effect [66]. Removal of 
neutrons from the 0/ 7 / 2 state, on the other hand, creates 
some neutron deficiency towards the interior, rather than 
the exterior, thus delaying the onset of the proton-skin 
effect in the dipole response. The above asymmetry be¬ 
tween neutron excess and neutron deficiency in the case 
of the Ni isotopic chain is obvious in the transition den¬ 
sities of the IS-LED, shown in Fig. 4. 

Finally, the splitting of the low-energy part of the re¬ 
sponse, from about 62 Ni to the shell closure, and the cor¬ 
responding unsmooth behavior of the IS-LED properties, 
constitutes a departure from what we observed, with the 
same model, in the Ca and Sn isotopic chains. In Ca, the 
filling of the vOf 7 / 2 subshell from 40 Ca to 48 Ca results 
in only a moderate and smooth change in the properties 
of the IS-LED. All involved nuclei are stable. In Sn, the 
filling of the j/2s — Id —0.g 7 / 2 —0/in/ 2 shell gives the same 
result, even though many of the nuclei are unstable. 

All the above observations underscore the relevance of 
shell structure in determining the El strength distribu¬ 
tion at low energies. Furthermore, shell effects such as the 
above have been found to influence the degree of correla¬ 
tion between the neutron-skin thickness and the amount 
of low-energy El strength [26, 64, 67, 68]. 


B. Pygmy states, polarizability, and asy-stiffness 

Next, we comment on the possible role of the symme¬ 
try energy and its slope L (asy-stiffness) in determining 
the low-energy El strength distribution, and in particu¬ 
lar how our results fit into the ongoing debate. A more 
profound discussion would require further quantitative 
analyses and, as we argue eventually, more accurate and 
complete data on 68 Ni. 

Our study of the Ca isotopic chain [17, 18] showed than 
no skin mode develops in 48 Ca at low energies (we cannot 
exclude one above threshold). The Gogny interaction, 
being asy-soft, predicts correctly that the IS-LED mode 
survives as such up to 48 Ca. Now we have found that the 
QRPA+D1S model underestimates the strength of the 
pygmy resonance, as measured around 10 MeV in 68 Ni, 
perhaps by a large amount. This could suggest that the 
functional is too soft - hardly a novel conclusion [61]. The 
same will be true for the SLy4 functional (also asy-soft), 
if the IS-LED and the neutron mode are bound. 
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Exp D1S SLy4 

208 Pb 
120 Sn 
68 Ni 

20.1(6) 18.78 20.075 
8.93(36) 8.45 9.40 

3.40(23) 3 3.79 4.03 


TABLE I. Dipole polarizability of three nuclei, in units of fm 3 , 
calculated within the QRPA-D1S and CRPA-SLy4 models, 
and measured values (“Exp”) from Refs. [7, 8, 38]. 


The effect of stiffness on the low-energy strength in 
208 Pb, 132 Sn, as well as 68 Ni, can be observed in the re¬ 
sults of Ref. [10]: In view of the measured strengths, the 
SkI3 functional (L = 100.5MeV) is too stiff, while the 
SLy4 (L = 46.0MeV) and SGII (L = 37.7MeY) func¬ 
tionals, and even more so the D1S (L = 22.4MeV), are 
too soft, in line with the conclusions of Ref. [61]. We 
note also that relativistic functionals, which are typically 
more asy-stiff [69, 70], predict strong neutron-skin modes, 
in the form of neutron-skin oscillations, for moderately 
asymmetric nuclei, thus overestimating their low-lying 
strength. 

We now point out that the softness effect is difficult 
to disentangle from shell effects. Extracting the asy- 
softness from data on a single nucleus, like 68 Ni, would be 
questionable, even though the correlations within a given 
model are found strong [5]. Rather, a different strategy 
to connect data with the symmetry energy and softness, 
and subsequent improvement of functionals, suggests it¬ 
self, akin to the spirit of Ref. [61]: The functional must be 
soft enough to reproduce the small El strength observed 
in the stable Ca and Sn isotopes below 8 — 9 MeV, as well 
as 208 Pb. At the same time, it should generate the correct 
amount of El strength in 68 Ni above threshold. Accu¬ 
rate data on 68 Ni are highly desirable. Furthermore, data 
below threshold would clarify the origin of the strength 
observed above threshold, with the help of our present 
results. 

An important related observation is that, although 
both models used here are rather soft, they overestimate 
the polarizability of 68 Ni. The same does not hold for 
other nuclei. In Table I we compare results obtained with 
the two models and measured values of ap in 68 Ni [38], 
120 Sn [8], and 208 Pb [7]: The QRPA+D1S model under¬ 
predicts the polarizability of 120 Sn and 208 Pb by 5% and 
7%, respectively, while the CRPA+SLy4 model describes 
accurately the polarizability of 208 Pb and overpredicts 
that of 120 Sn by 5%. Both models overestimate the po¬ 
larizability of 68 Ni by 10% or more. If we trust the model 
systenratics, namely that such soft models cannot over¬ 
estimate the dipole polarizability by such a factor, then 
the discrepancy represents a large amount of El strength 
potentially missing from the data. 3 


3 It is reported [47] that a revised measured value of 3.88(31) fm 3 


We iterate that, besides symmetry and isovector prop¬ 
erties of the nuclear energy-density functional, additional 
factors determine the strength distribution, which in¬ 
clude shell effects and phonon coupling, as well as the 
nucleon effective mass in the case of mean-field function¬ 
als and the convergence of the calculations just above 
particle threshold. We conclude, in any case, that more 
accurate and complete data are needed for 68 Ni, than 
reported in Ref. [38], before the polarizability of this nu¬ 
cleus can be used to calibrate theoretical models or con¬ 
strain the equation of state with confidence. 

V. SUMMARY AND CONCLUSIONS 

We have calculated the dipole response of even 
Ni isotopes using the QRPA+D1S Gogny model and 
the CRPA-|-SLy4 Skyrrne model, which predict similar 
nuclear-matter properties of relevance (symmetry energy 
and softness, nucleon effective mass). We have analyzed 
thoroughly the dipole transitions at low energy. Our 
main conclusions can be summarized as follows: 

• The two models generally agree in their predictions 
regarding the various transitions below 11 MeV, ex¬ 
cept for additional continuum strength in the case 
of CRPA. Above 11 MeV and in the region of the 
GDR there are discrepancies in the IV channel, at¬ 
tributable to the different TRK enhancement fac¬ 
tor. 

• The summed El strength of stable isotopes at low 
energies is well reproduced. 

• In all isotopes, a strong IS low-energy dipole mode 
(IS-LED) is obtained around 10 MeV with the 
QRPA+D1S model. Given the systenratics of 
our model, we expect it to actually occur around 
7 MeV. In the stable and moderately neutron-rich 
isotopes it is expected to be a bound transition. 

• We showed that shell structure plays an important 
role in determining the low-energy El strength dis¬ 
tribution. Furthermore, we conjectured that the 
isospin structure of low-energy strength, observed 
in stable nuclei, is determined by the difference be¬ 
tween the IS-LED energy and the 1 Hu> value. 

• In 68 Ni and neighboring isotopes, besides the IS- 
LED, one more coherent IS transition is found, ly¬ 
ing a couple of MeV above the IS-LED. Complete 
data are needed to determine whether there is a bi- 
rnodal structure below threshold, or whether we can 
associate this second peak with the one detected 
above threshold [38, 48]. 


is obtained by extrapolating a Lorenzian fit to the data [38] be¬ 
yond the measured energy window. That could lift some of the 
discrepancy observed here. 
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• Beyond the shell closure, 78 Ni, the neutron-rich iso¬ 
topes are found highly polarizable, due to valence- 
neutron transitions. 

• We showed that the relation between the polar¬ 
izability and the skin thickness along an isotopic 
chain is not straightforward, because they are in¬ 
fluenced by shell effects in different ways or degrees. 

• The QRPA+D1S model underestimates the thresh¬ 
old El strength in 68 Ni. Furthermore, if both 
inodes identified above (IS-LED and neutron mode) 
are bound, the model fails to explain the peak ob¬ 
served above threshold. One reason for these prob¬ 
lems must be the asy-softness of the functional, or 
the residual interaction in the IV channel, pushing 
too much strength into the GDR. Other reasons 
may include the nucleon effective mass, effects be¬ 
yond QRPA, especially phonon coupling [21 24], 
and the continuum treatment, the latter of which 
we examined. 

• We found that a converged treatment of continuum 
threshold states, as in the CRPA model, leads to 
additional threshold strength in neutron-rich nu¬ 
clei. In 78 Ni its magnitude could be significant. 

• Although the present functionals are soft, mostly 
underestimate the dipole polarizability of stable nu¬ 
clei and underestimate the El strength in 68 Ni, 
they overestimate the polarizability of 68 Ni. This 
leads us to infer that strength is missing from the 
68 Ni data reported in [38]. 

• We conjectured, pending further systematic stud¬ 
ies, that accurate and complete data on 68 Ni not 
only could clarify the above issues, but provide a 
complementary constraint to the asymmetry stiff¬ 


ness from below, while data on stable nuclei could 
constrain it from above. 

In the future, it will be very interesting to confirm ex¬ 
perimentally the dramatic role of shell closure revealed 
in this and other theoretical studies. This would require 
data along isotopic chains, in particular nuclei on either 
side of a shell closure, and in the same energy range, 
especially in relation to the particle emission threshold. 

In short, the present study of the Ni chain, with its par¬ 
ticular shell structure and many isotopes, has helped us 
identify new phenomena and factors which determine the 
low-energy dipole strength. New data above and below 
particle-emmision threshold should enable us to test our 
conclusions and provide important feedback for nuclear 
structure theory. 
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